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INTRODUCTION 

i ■ acolication of CLS and NESSUS 

This paper describe t o P t he probabilistic structural 

family of 5°®P^ e | SME g paC e propulsion system components, 
analysis of four SSME space pr P enviro J ments that are 

- «- 

£££r“S..? S n retirement 

Srm^irl^^^ee^^lXy used in the past in 

deterministic analysis studies. 

PROBABILISTIC TURBINE BLADE STATIC AND MODAL ANALYSIS 

A high P r ® 3 ^ r %^i e t if°notaroi n Sinteen ^Ldom 31 ^^ 
variables^were ^considered in the analysis J° Win9 

p r :^tirn 9 :i?Jarior^ni?r:rrhrtna?ysirpne f ctin g 

Inaction . KI^PpWcatl^LIn^^ 

probabilistic static «alysis as can 

heini e de>i9™"9 1 "° r ® ndurable components. (Reference 1). 

sssi’s 

"r”.r:s*:p'.L ss~ *“■ 

quantitative probability estimates of interference »rtn 
engine orders. The first ten variables listed in Table 
were used in this probabilistic analysis. 

The results from a mean value first order and the mor ® 
accurate advanced mean value first order results are shown 
In I^ule 4 for mode 1. It points out *he asymmetrical 
nature of the distribution and the capability of the 
probabilistic tools in predicting it. The computed 
coefficient of variation of the natural f: requ< encies is 
consistant with the actual experience. If the freque y 
range is a criteria for acceptance or rejection °f 
thil analysis technics can then be used to calculate 
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target design point 


that will reduce the rejection rate 


(Reference 2) . 

VIBRATION ANALYSIS OF DUCTS SUBJECTED TO UNCERTAIN LOADS 


The high pressure oxidizer turbopump discharge duct is used 
in this study (Figure 5) . The duct has three attach points 
subjected to both random and sinusoidal vibration, a 
multisupport excitation problem. The variation in the 
vibratory loads is attributed to engine system duty cycle 
operation, engine system hardware variations, and local 
component variations within the turbopump or combustors. A 
total of thirty eight random variables were used in this 
analysis (Sbl e 2) A typical result obtained in the form of 
a cumulative distribution function at a typical node for 
bending moment in the y direction is shown in figure 6. The 
analysis points to a way of designing structures subjected 
to a large number of excitation sources without exercizing 
undue conservatism (Reference 3) . 

PROBABILISTIC MATERIAL NONLINEAR ANALYSIS OF LOX POST 


A main injector element of SSME is used in this ®tv» d y . 
(Figure 7). The response variable of interest is the cyclic 
strlin range (including the elastic and plastic portions of 
the stress strain curve) . This response quantity is one of 
the primary drivers in determining Low Cycle Fatigue life o , 
the component. The dominant loading on the component is the 
differential wall temperature across the Lox Post wall of 
approximately 1000 degrees rankine. It was the only loading 
considered in this analysis. 


The random variables that affect the local temperature f ^ eld 
are shown in Table 3. The probabilistic analysis involved a 
full two duty cycle incremental nonlinear analysis and 
constructiing a response surface linking effective strain 
range to random variables . The probabilistic analysis 
results based on the response surfaces at several locations 
are shown in Table 4 and the corresponding sensitivity 
factors in Table 5. The analysis demonstrates a methodology 
of linking global system variables to local response 
variables that can eventually be extended to calculate 
damage, all in a probabilistic domain (Reference 4) . 


PROBABILISTIC BUCKLING LOAD ANALYSIS 


The structural liner in the SSME two duct hot gas manifold 
design was the subject of this study. The shell is doubly 
curved with five distinct zones of thicknesses (Figure 8) . 
One of the primary design requirement for the liner is to 
have adequate margin against buckling failure. Unlike the 
previous examples cited above this application deals with 
estimation of the strength variable. As a first step in this 
process a probabilistic linear buckling analysis was 
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conducted. A more rigorous material and geometric 
probabilistic collapse load analysis would be more accurate 
but the computational effort will also be significantly 
larger . 

In this study the thickness in the five zones of the shell 
were considered as independent random variables but within 
each zone, the thickness variation was considered to be 
fully correlated. The probabilistic analysis method used is 
a response surface approach. The resulting computed 
cumulative distribution function is shown in Figure 9. 

Thus knowing the distribution of the buckling strength and 
the distribution of differential pressure the probability of 
buckling failure based on the linear eigen value analysis 
can easily be calculated. 

SUMMARY 

Application of the probabilistic analysis tools developed in 
CLS AND PSAM contracts to a select SSME components has been 
successfully demonstrated. The scope and size of the 
application prove the viability and usefulness of the tools 
and methods developed to practical design in terms of 
designing more durable components. 
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TABLE I 

LIST OF RANDOM VARIABLES USED IN 
TURBINE BLADE ANALYSIS 


_ RANDOM VA RIABLE 
NO DESCRIPTION 

1 MATERIAL AXIS 
ABOUT Z 

2 MATERIAL AXIS 
ABOUT Y 

3 MATERIAL AXIS 

ABOUT X 

4 ELASTIC MOOULUS 


5 POISSON'S RATIO 

6 SHEAR MOOULUS 

7 MASS OENSITY 

8 GEOMETRIC LEAN 
ANGLE ABOUT X 

9 GEOMETRIC TILT 
ANGLE ABOUT Y 

10 GFOMETRIC TWIST 
ANGLE ABOUT Z 

11 MIXTURE RATIO 
LIQUID HYDROGEN/ 
LIQUID OXYGEN 

12 FUEL INLEI 
PRESSURE 

13 OXIOIZER INLET 
PRESSURE 

14 FUEL INLET 
TEMPERATURE 

15 OXIDIZER INLET 
TEMPERATURE 

16 HPFP EFFICIENCY 

17 HPFP HEAD 

COEFFICIENT 

18 COOLANT SEAL 
LEAKAGE FACTOR 

19 HOT GAS SEAL 
LEAKAGE FACTOR 


MATERIAL 


VAR1A1 IONS 


ELASTIC 

PROPERTY 

VARIATIONS 

MASS 

VARIATIONS 


GEOMETRY 

VARIATIONS 


INDEPENDENT 


VARIATIONS 
DEPENDENT 
LOADS ARE 
TURBINE BLADE 
PRESSURE 
TEMPERATURE 


AND SPEED 


LOCAL 


6E0HETRY 


FACTORS 


( FEM 
(QUANTITIES 
I AFFECTED 

I 

| MATERIAL 

I 

| ORIENTATION 

I 

| ANGLES 


ELASTIC 

CONSTANTS 


NOOAL 

COORDINATES 


PRESSURE 


TEMPERATURE 


CENTRIFUGAL 


TEMPERATURE 


-0.087266 
RAO IANS 

-0.034907 

RADIANS 

+0.052360 
RAOIANS 
18.3BE6 PSI 
(126.22E9 Pa) 

0.386 

18.63E6 PSI 
(12B.45E9 Pa) 

0.805E-3 


30.00 PSI 
(2.068E5 Pa) 

100.00 PSI 
(6.894E5 Pa) 

38.5* R 
(21.39* K) 

167.0* R 
(92.78* K) 


1.0237 


STANOARO 

DEVIATION 

0.067544 

RAOIANS 

0.067544 

RADIANS 

0.067544 
RADIANS 
0.4595E6 PSI 
(3.168E9 Pa) 

0.00965 

0.46575E6 PSI 
13.223E9 Pal 
0.493E-5 


0.14 

DEGREES 

0.14 

DEGREES 

0.30 

DEGREES 


5.00 PSI 
( .344E5 Pa) 

26.00 PSI 
(1.793E5 Pa) 

0.5* R 
(0.278* K) 

1.33* R 
(0.739* K) 

0.008 


0.008 
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Location A 



Figure 1. Turbine Blade Finite Element Model 
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Figure 2. CDF For Effective Stress at Location A 
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Figure 3. CDF For Effective Stress at Location B 


UOOE 1 - TMtGENlM. HEX MODE 



FREQUENCY, HERTZ 
□ MYFO + ADiMTO 


Figure 4. CDF For Mode 1 Frequency Using HVFO and ADMVFO 
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Table 2. HPOTP Discharge Duct Analysis 
Random Variable Statistics 


Sequence 

No. 

Random Variable Description 

Mean 

CV 

Dist. Type 

1) 

Zone G - X axis, PSD power 

222.0 

0.73 

Log Normal 

2) 

Zone G - Y axis, PSD 

73.5 

0.808 

Log Normal 

3) 

Zone G - Z axis, PSD 

73.5 

0.808 

Log Normal 

4) 

Zone A - X axis, PSD 

22.5 

0.20 

Log Normal 

5) 

Zone A - Y axis, PSD 

54.0 

0.20 

Log Normal 

6) 

Zone A - Z axis, PSD 

69.5 

0.2 

Log Normal 

7) 

Oxidizer Pump Speed 

2940.53 

0.014 

Log Normal 

8) 

Fuel Pump Speed 

3707.08 

0.01 

Log Normal 

9) 

Damping 

0.033 

0.15 

Normal 

10) 

Zone A Main Injector 
Oxidizer Pump Sine Amplitudes 

X direction IN 

0.30 

0.4 

Log Normal 

11) 

2N 

0.30 

0.15 

Log Normal 

12) 

4N 

1.5 

0.3 

Log Normal 

13) 

Y direction IN 

0.60 

0.5 

Log Normal 

14) 

2N 

0.70 

0.40 

Log Normal 

15) 

4N 

2.6 

0.3 

Log Normal 

16) 

Z direction IN 

0.5 

0.45 

Log Normal 

17) 

2N 

0.70 

0.20 

Log Normal 

18) 

4N 

0.70 

0.20 

Log Normal 

19) 

Fuel Pump Sine Amplitudes 
X Direction IN 

0.35 

0.3 

Log Normal 

20) 

Y Direction IN 

0.80 

0.35 

Log Normal 

21) 

Z Direction IN 

1.20 

0.3 

Log Normal 

22) 

Zone G - Oxidizer Turbopump 
Oxidizer Pump Sine Amplitudes 
X Direction IN 

1.35 

1.0 

Log Normal 

23) 

2N 

1.50 

0.5 

Log Normal 

24) 

3N 

1.10 

0.45 

Log Normal 

25) 

4N 

11.0 

0.25 

Log Normal 

26, 27) 

Y&Z Direction IN 

1.9 

0.9 

Log Normal 

28. 29) 

2N 

1.6 

0.6 

Log Normal 

30, 31) 

3N 

0.75 

0.3 

Log Normal 

32, 33) 

4N 

5.5 

0.6 

Log Normal 

34) 

Fuel Pump Sine Amplitudes 
X Direction IN 

0.65 

0.35 

Log Normal 

35, 36) 

Y&Z Direction IN 

0.45 

0.4 

Log Normal 

37, 38) 

Y&Z Direction 2N 

0.45 

0.4 

Log Normal 

Note: 1) Power units are in g 2 

2) Pump speed units are in radians/second 

3) Sinusoidal amplitude units are in g 
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FIGURE 5. DUCT MODEL WITH ELBOWS, VALVES, 
ATTACHMENTS AND SECONDARY LINES. 



FIGURE 6. CDF OF A RESPONSE BASED ON 
MVFO AND AMVFO METHODS. 
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FIGURE 7. LOX POST AND ITS AXISYMMETRIC FINITE ELEMENT MODEL. 


TABLE 3 

RANDOM VARIABLES USED IN LOX-POST ANALYSIS 


Random Variable 

Mean 

Standard 

Deviation 

Distribution 

Material yield stress (ksi) 

175.0 

8.75 

Normal 

Hot-gas temperature (R) 

1654.70 

26.6407 

Normal 

Coolant temperature (R) 

191.643 

4.21019 

Normal 

Hot-gas flowrate (lbm/sec) 

167.249 

1.0928 

Normal 

Cooiant flowrate (lbm/sec) j 

929.918 

4.31211 

Normal 

Mixture ratio 

0.948012 

0.0184211 

Normal 

Heat-shield-gap factor 

0.47 

0.235 

Lognormal 

Hot-gas film coefficient 

1.0 

0.1 

Normal 

Coolant film coefficient 

1.0 

0.08 

Normal 
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Node 

Median 

Mean 

Standard 

Deviation 

56 

0.004143 

0.004821 

0.0002613 

78 

0.006759 

0.008360 

0.0002275 

84 

0.002114 

0.004175 

0.0001661 

117 

0.002990 

0.004757 

0.0001454 

122 

0.001186 

0.003831 

0.0001473 

175 

0.002945 

i 0.00513f 

0.0001430 

181 

0.00229C 

j 0.002401 

O 

o 

o 

o 

o 

ON 

4^ 

o 


°f 

V ariation 
0.054 
0.027 
0.040 
0.031 
0.038 
0.028 
0.027 


Random Variable 
Hot-gas temperature 
Coolant temperature 
Hot-gas flowrate 
Coolant flowrate 
Mixture ratio 
Shield-gap factor 
Hot-gas film coefficient 
Coolant film coefficient 


Node 56 
0.455 
0.014 
0.062 
0.000 
0.022 
0.793 
0.399 
0.002 


Node 78 
0.796 
0.026 
0.075 
0.003 
0.034 
0.115 
0.587 
0.031 


Node 84 
0.771 
0.024 
0.111 
0.002 
0.034 
0.071 
0.620 
0.042 


Node 117 
0.802 
0.023 
0.070 
0.002 
0.035 
0.028 
0.590 
0.024 


Node 122 
0.793 
0.023 
0.074 
0.002 
0.034 
0.003 
0.603 
0.026 


Node 175 
0.860 
0.019 
0.045 
0.001 
0.028 
0.011 
0.507 
0.005 


0.786 

0.022 

0.108 

0.003 

0.031 

0.017 

0.605 

0.052 


Range between 0 and 1. Larger values indicate a greater influence of the random variable on 
the response. 
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FIGURE 8. FINITE ELEMENT MODEL OF THE 
STRUCTURAL LINER. 


CUMULATIVE DISTRIBUTION FUNCTION 



PRESSURE PSI 

O UVFO + MONTEORLO 

FIGURE 9. CDF OF THE FIRST BUCKLING MODE PRESSURE. 
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